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We report a high precision search for orbital-like magnetic order in the pseudogap region of 
I^-^Sr^CuCU single crystals using zero-field muon spin relaxation (ZF-^iSR). In contrast to pre- 
vious studies of this kind, the effects of the dipolar and quadrupolar interactions of the muon with 
nearby nuclei are calculated. ZF-piSR spectra with a high number of counts were also recorded to 
determine whether a magnetically ordered phase exists in dilute regions of the sample. Despite these 
efforts, we find no evidence for static magnetic order of any kind in the pseudogap region above the 
hole-doping concentration p = 0.13. 

PACS numbers: 74.72.Gh, 75.25.Dk, 76.75.+i 



A distinctive feature of high transition tempera- 
ture (T c ) cuprate superconductors is the pseudogap re- 
gion that exists above T c and over a wide range of 
doping. For some time there has been much de- 
bate on whether the pseudogap is a manifestation of 
a phase transition. Recently it has been demon- 
strated that the onset of the pseudogap in optimally- 
doped Pbo.55Bii.5Sri.6Lao.4Cu06+5 at a temperature 
T* w 3.5T C is likely a sign of a phase transi- 
tion to a non-superconducting state with a broken 
symmetry^ Although this is not necessarily a mag- 
netically ordered state, proposed orders for the pseu- 
dogap state include time-reversal symmetry breaking 
phases that contain ordered circulating orbital currents, 
which either break 2 ^ 3 - or preserve^ translational sym- 
metry. The strongest experimental evidence for an 
orbital-current phase are the the observations of an un- 
usual translational-symmetry preserving magnetic or- 
der in YBa2Cu30 y and HgBa2Cu04+5 by spin-polarized 
neutron diffraction, which bear some resemblance to 
the ordered 0/j circulating-current phase proposed in 
Ref. 0- It is worth pointing out, however, that by ex- 
trapolation the onset of this orbital-like magnetic order 
is expected to occur near T c at optimal doping, and hence 
its relationship to the phase transition reported in Ref. 1 
is unclear. Furthermore, to date such orbital-like mag- 
netic order has not been observed in any cuprate beyond 
a hole doping of p = 0.135, whereas the 0/j phase is 
predicted to persist up to p~0.19. 

In contrast to the experimental techniques used in the 
above studies, local probes of magnetism, such as nu- 
clear magnetic resonance (NMR) and zero-field muon 
spin relaxation (ZF-/^SR) provide information on the 
magnetic volume fraction. Unfortunately, 89 Y NMR 
experiments on Y2Ba4Cu70i5_5f U ^ Zeeman perturbed 
nuclear quadrupole resonance (NQR) measurements of 
YBa 2 Cu40sri and ZF-/xSR experiments on YBa2Cu30 y 
(Ref.QlGJ) and La 2 _ K Sr 2; Cu04 (Ref.|]J) have found no 
evidence for the onset of magnetic order at the pseu- 
dogap temperature T*. One exception is the finding 



of anomalous magnetic order by ZF-yLtSRji 3 - in the same 
YBa 2 Cu306.6 single crystal studied in Ref. H. The mag- 
netic order is characterized by an onset temperature and 
an average local dipolar magnetic field that are in quan- 
titative agreement with the orbital-like magnetic order 
detected by polarized neutron diffraction. Yet the ZF- 
/iSR measurements clearly show this form of magnetic 
order existing in only about 3% of the sample, suggest- 
ing that it is associated with a minority phase in lower 
quality samples. 



Orbital-like magnetic order has also been observed in 
Lai.9i5Sro.o85Cu04 by polarized neutron diffraction, 15 
but is less pronounced than the long-range magnetic or- 
der that has been reported in the pseudogap regions of 
YBa 2 Cu 3 O y (Ref.i,i) and HgBa 2 Cu0 4+ 5 (Ref.0). In 
particular, the magnetic order in Lai.gi5Sro.o85Cu04 is 
short range, two (rather than three) dimensional, and oc- 
curs at a temperature far below T* . Yet no such magnetic 
order was observed in a ZF-//SR study of x > 0.13 samples 
by MacDougall et alr^ Nevertheless, the subtle nature of 
the magnetic order observed in Lai.gi5Sro.o85Cu04 war- 
rants a more precise ZF-fj.SH search to verify its existence 
and determine whether it is an intrinsic property of the 
pseudogap phase of La2_ a; Sr a; Cu04. 



Here we report two significant advances in the applica- 
tion of ZF-/iSR to search for magnetic order in the pseu- 
dogap region of La2- a: Sr :l: Cu04. First, we accurately de- 
termine the interactions of the positive muon (/i + ) with 
the nuclear spin system, allowing us to identify any resid- 
ual relaxation of the ZF-//SR spectrum that could be as- 
cribed to static magnetic order. In doing so we have 
accurately identified the muon stopping site. Second, we 
have acquired ZF-/1SR spectra of higher statistics than 
in previous works, enabling a search for dilute or short- 
range magnetic order. 
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I. EXPERIMENTAL DETAILS 



The experiments were performed on platelet-like single 
crystals of I^-zSr-j-CuO^ cut from a travelling-solvent 
floating zone (TSFZ) growth rod. The TSFZ growth pro- 
cedure that was followed is decribed elsewhere*^ The 
single crystals cut from the TSFZ rod were annealed at 
800 to 900 °C in an oxygen partial pressure to remove 
oxygen defects in accordance with the oxygen nonstoi- 
chiometry of La2_ ;c Sr 2 ,Cu04jiZ Strontium concentrations 
greater than x = 0.125 were chosen to ensure the absence 
of static antiferromagnetism or the spin-glass-like mag- 
netism previously observed by ZF-/iSR in lower doped 
samples! 18 ' 19 Magnetic susceptibility measurements of 
the bulk superconducting transition temperature by a 
superconducting quantum interference device (SQUID) 
yield T c = 37.6, 37.3, 28 and 17 K for the x = 0.15, 0.166, 
0.216, and 0.24 samples, respectively. 

The ZF-/iSR measurements of the La2- K Sr :E Cu04 sin- 
gle crystals were performed on the M15 and M20B surface 
muon beam lines at TRIUMF. Positive muons implanted 
into the sample Larmor precess about the local field B 
and decay according to fi + — > e + + v e + 9^, with a mean 
life time r ~ 2.2 fxs. The ZF-/xSR signal is generated 
from detection of the decay positrons, which are prefer- 
entially emitted along the muon spin direction. The sam- 
ples were mounted with the crystallographic c-axis of the 
La2- 2; Sr a ;Cu04 single crystals parallel to the muon beam 
momentum. The initial muon spin polarization P(0) was 
oriented perpendicular to the c-axis using a Wien filter. 
This has the advantage that neither positron detector di- 
rectly faces the incoming muon beam. In this geometry 
the ZF-//SR "asymmetry" spectrum is defined as the dif- 
ference between the number of decay positrons sensed by 
scintillator detectors positioned above (A) and below (B) 
the sample, divided by the sum of the counts in these two 
detectors 



N A -N B 
Na + Nb 



A(t) = aP(t) . 



(1) 



where a < 1/3 is the initial asymmetry (dependent on the 
energy of the decay positrons and several experimental 
factors) and P(t) is the time evolution of the muon spin 
polarization. The latter is modeled by an appropriate 
relaxation function G(t) 



Pit) = G(t) cos^Bt) , 



(2) 



where 7^ is the muon gyromagnetic ratio and B is the 
average internal magnetic field sensed by the muon. Note 
that B = in the absence of magnetic order. 



II. EXPERIMENTAL RESULTS 

In previous studies of I^-zSr^CuO^ the contribution 
of the nuclear moments to the ZF-/1SR signal was as- 
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FIG. 1: (Color online) ZF-/1SR signal of single crystal 
Lai.784Sro.2i6Cu04 at T = 30 K recorded with the initial muon 
spin polarization P(0) perpendicular to the c-axis. The solid 
green curve is a fit of the data below t = 6 fis to Eq. <(3j mul- 
tiplied by an asymmetry factor a. The solid black curve is a 
fit to Eq. ©. 



sumed to be described by a static Gaussian Kubo-Toyabe 
(KT) relaxation functio n 20 ! 21 



G KT W = i + ^(l-A 2 i 2 )exp 



--A 2 t 2 
2 



(3) 



where A is the second moment of the local magnetic field 
distribution at the muon site. The relaxation function 
Gkt(^) assumes that the muon is immobile, the local 
magnetic fields acting on the muon spin are static, and 
that the local magnetic field distribution is isotropic and 
Gaussian. The nuclear dipolar fields sensed by the posi- 
tive muon /i + are usually static, because correlation times 
of the nuclear moments are generally much longer than 
the muon life time. Assuming that the nuclear dipolar 
fields acting on the muon spin are randomly oriented and 
that each muon sees a unique local field over the duration 
of its life time, 1/3 of the muon spins will be parallel to 
the field and not evolve in time. Consequently, the re- 
covery of the muon spin polarization to 1/3 of its initial 
value at late times is characteristic of static fields. How- 
ever, in a single crystal the contribution of the surround- 
ing nuclear moments can substantially deviate from the 
random field approximation, with no recovery of P(t) to 
1/3. For example, Fig. [TJ shows that the static Gaussian 
Kubo-Toyabe function does not fully describe the ZF- 
/iSR spectrum of single crystal Lai.784Sro.2i6Cu04 for 
the case where P(0) is perpendicular to the c-axis — de- 
spite being in a region of the phase diagram where only 
the nuclear dipole moments are expected to contribute to 
the ZF relaxation function. Of particular note, the muon 
spin polarization dips below zero, whereas Gkt(0 does 
not. To circumvent this problem the ZF-fiSK signals for 
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FIG. 2: (Color online) (a) Comparison of the ZF-/1SR signals 
of Lai.834Sro.i66Cu04 at T = 2.5 K and T=158 K measured 
over a 10 /is time range with P(0) perpendicular to the c-axis. 
(b) The early time ZF-/1SR signal of Lai.834Sro.i66Cu04 at 
T — 2.5 K shown over the first 1 /is. In both figures the solid 
curves are fits to the relaxation function of Eq. ([4]) . 



FIG. 3: (Color online) Results of fits to Eq. 0}. (a) Temper- 
ature dependence of B for all four samples, (b) Dependence 
of the relaxation rate A (solid squares) and the exponent K 
(open circles) on Sr content x. Note, both A and K are tem- 
perature independent fit parameters. 



La2- 2; Sr 2 ;Cu04 in Ref. were truncated at t = 8 /is to 
facilitate fits to Eq. (0. 

To more accurately portray changes in the functional 
form of the ZF-/1SR signal as a function of temperature 
and Sr content x {i.e. hole doping concentration), here 
we fit the ZF-/tSR spectra to a simple phenomenological 
power-exponential relaxation function, such that 

P(t) = cxp [-(At)*] cos[ 7m S(T)<] , (4) 

with the relaxation rate A and the exponent K treated as 
temperature independent quantities, and the average lo- 
cal magnetic field B being the only fit parameter allowed 
to vary with temperature. Figure Ufa) shows representa- 
tive time spectra for Lai.834Sro.i66Cu04 that are well de- 
scribed by Eq. © over the entire 10 /is time range. A fit 
to Eq. (0 is also shown for Lai.784Sr .2i6CuO4 in Fig. Q] 
The evolution of the ZF-/1SR signal with temperature is 
reflected in the temperature dependence of B shown in 
Fig. 0a). Consistent with the findings of MacDougall et 
al.r^ the ZF-/1SR signal does not exhibit a temperature 



dependence characteristic of a magnetic phase transition 
in any of the samples. At T = 200 K, the value of B 
is clearly reduced, but by the same amount at x = 0.15 
and £ = 0.24. Hence, the reduction of B at high temper- 
atures is likely caused by muon diffusion, whereby the 
mobile /t + sees a time-averaged field over its life time. 
This same conclusion was reached in Ref. 

The dependence of A and K on the Sr content x is 
shown in Fig. 0b). There is some increase in the values of 
both parameters with increasing x, but a clear reduction 
of A and K for the 2: = 0.24 sample. The former behav- 
ior may be the result of antiferromagnetic fluctuations, 
dilute regions of static magnetism and/or the relaxation 
caused by the Sr nuclei. While the precise source is un- 
clear, there is no evidence for static magnetic order in 
any of these samples. The smaller values of A and K at 
a; = 0.24 indicate another contribution to the ZF-/iSR sig- 
nal. This is likely paramagnetic moments that are known 
to be present in heavily-overdoped La2- a: Sr : j : Cu04 above 
x^0.19i^2r— The onset of the Curie-like paramagnetism 
is the probable cause of the slightly larger values of B at 
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x = 0.216 and x = 0.24 in Fig.^a). 

In Ref. [l3l three distinct magnetic components were 
observed in the ZF-/jSR signal of a large YBa2Cu306.6 
single crystal. Neutron scattering experiments on this 
same single cystal show that the sample contains mag- 
netic order associated with the impurity "green phase" 
Y 2 BaCu05r& and two additional kinds of unusual mag- 
netic orderi&2& While one of the latter two magnetic 
components is manifested as a slow relaxing component 
that is also observed in higher quality samples, the other 
two forms of magnetic order are discernible as small- 
amplitude, rapidly-damped oscillatory components in the 
early time range of the ZF-/iSR signal. The small ampli- 
tudes indicate that the magnetic orders are confined to 
small volume fractions of the YBa2Cu30g.6 single crys- 
tal. However, Fig. HJb) shows that there are no such 
oscillatory components in the early time ZF-/uSR signal 
of Lai.834Sro.i66Cu04. We have also examined the early 
time spectra of the other samples at various temperatures 
above T = 2.3 K, and likewise find no evidence for any 
kind of short-range or dilute magnetic order. 

Next we calculate the contribution of the nuclei to the 
polarization function P{t), in an effort to fully account 
for the observed ZF-/iSR spectrum. 



III. CALCULATION OF NUCLEAR-INDUCED 
RELAXATION 

Here we describe a general numerical method for calcu- 
lating the ZF relaxation function resulting from the dipo- 
lar magnetic and quadrupolar electrostatic interactions of 
the n + with an arbitrary number of neighboring nuclei. 
We consider only the interactions of the nuclei with the 
muon, and ignore interactions amongst the nuclei them- 
selves. This approximation is justified by the size of 7^, 
which is about an order of magnitude larger than the 
gyromagnetic ratios of the nuclei in La2- :c Sr. E Cu04. We 
also assume that the electric field gradient (EFG) at each 
nuclear site is due to the Coulomb field of the unscreened 
positively charged muon, and the non-symmetric charge 
distribution of the crystal itself (i.e. the crystal EFG). 

The interactions between the muon spin (S — 1/2) 
and N surrounding nuclei of spin /, as well as the ef- 
fect of the crystal EFG, is described by the following 
Hamiltonia n 27 ! 28 
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H 



H, 



(5) 



where 
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In the above equations, the term Hj 3 is the dipole- 
dipole interaction between the positive muon and the j th 
nucleus, is associated with the quadrupolar energy of 
the nuclear spin Ij due to the EFG generated by the pos- 
itive muon, iij is the unit vector pointing in the direction 
along the straight line that connects the muon to the j th 
nucleus located a distance fj away, and 7^ and jj are the 
gyromagnetic ratios of the muon and nuclei, respectively. 
The quadrupolar coupling constant uj®j is proportional 

to l/r 3 ^ The term represents the quadrupolar en- 
ergy of the nuclear spin due to the crystal EFG, with 
a quadrupolar coupling constant u/^-. The constant rj is 
an asymmetry parameter which specifies the symmetry 
of the crystal EFG around the nucleus^ In our calcula- 
tion, the values of the quadrupolar coupling and the 
asymmetry parameter rj for La and Cu are taken from the 
literature. In particular, uj^- — 2iivq with vq = 34.0, 31.0 
and 6.40 Hz for 63 Cu, 65 Cu and 139 La, respectively^^ 
At Cu sites, 77 = 0.03, 32 while at La sites, 77 = 0.02^. In 
addition, the weighted averages of the two isotopes 63 Cu 
and 65 Cu are used for the gyromagnetic ratio 7c u an d 
the nuclear quadrupole moment Qcu of the Cu nuclei. 

Since the only stable isotope of Sr with nonzero spin 
is 87 Sr with a natural abundance of 7%, the Sr nuclei 
can be neglected in the calculations. We note that this is 
consistent with the study by MacDougall et al.^ 14 which 
showed that the relaxation rate of the ZF-/iSR signal of 
La2- a; Sr 2 ;Cu04 varies little in the range 0.13 < x < 0.30. 
Consequently, it is sufficient to perform our calculations 
for a positive muon residing in a single crystal of the 
parent compound La2Cu04. 

With knowledge of the Hamiltonian [Eq. ([5])], the time 
dependence of the muon spin polarization may be calcu- 
lated from the density matrix of the spin system as 

P(t) = Tr {p(0)e^ [a, ® (®f =1 l^)] e"^} . (9) 

Here oy, = cy, ■ P(0) = a x sin 9 cos + a y sin 6 sin j3 + 
a z cos 9 is the projection of the muon spin along the di- 
rection of the initial polarization P(0), with 9 and f3 the 
polar and azimuthal angles, respectively, between P(0) 
and the crystallographic axes (see Fig. [4}. Note that 
from Eq. (JSj> , the crystal EFG sets the quantization axis 
to be along the c-axis of La2-xSr a Cu04. The initial den- 
sity matrix is p(0), Dj = 21 j + 1 is the spin degeneracy for 
the j th nuclear spin, and l n is the nxn identity matrix. 

For a system in which the nuclear spins are randomly 
oriented, the initial density matrix is 



H% - tw% [(Ij • n 3 ) • (Ij ■ iij) - 1(1 + l)/3] , (7) 
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FIG. 4: (Color online) The muon D site (0.120, 0, 0.219) in 
the tetragonal unit cell of La2Cu04, which best describes the 
ZF-^iSR signal of La 2 _ ;c Sr :c Cu04 at x > 0.13 (see Fig. 0. The 
coordinate values are in multiples of their respective lattice 
constants a = 3.80 A, fe = 3.80 A and c= 13.12 A.— The small 
shaded rectangle visually emphasizes that the muon at the D 
site near the apical O resides on the a-c face of the unit cell. 
The lower left of the picture shows the 5 La nuclei and the 1 
Cu nucleus included in the calculation of P(t) for this muon 
site (Two of the La atoms are labeled 1 and 2 to identify them 
in the full crystal structure). 
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FIG. 5: (Color online) (a) Approximation method calcula- 
tions of the muon spin polarization function Pit) for different 
muon sites and with the initial muon spin polarization P(0) 
perpendicular to the c-axis and oriented at 45° with respect 
to the a-axis. The spatial coordinates of the muon sites and 
the nuclei included in each calculation are shown in Tabled] 
The sites Mi, Ti, T L > and Ui were previously considered in 
Ref. HH. (b) P(t) for the same muon site, but with different 
nuclei used in the calculation (see Table U]). 



where 



D 



N 
3=1 



(11) 



is the dimensionality of the Hamiltonian matrix. The 
time evolution of the muon spin polarization P(t) may 
be determined exactly by diagonalizing the DxD Hamil- 
tonian. The main limitation of this approach is that 
it becomes very computationally expensive as the num- 
ber of nuclei (and hence D) increases. In this work the 
largest systems studied include 10 nuclei, corresponding 
to D = 2 26 «6.71xl0 7 . We found that for these large val- 
ues of D, it was necessary to use an approximate method 
to determine P(t), and hence we used a method devel- 
oped by Celio to study P{t) in copper^ This approxi- 
mation method is based on the Trotter formula and uti- 
lizes the random phase approximation. We verified good 
agreement with the exact calculation for some of our cal- 
culations with less nuclei. 



The calculation of P(t) is strongly dependent on the 
orientation of P(0) with respect to the c-axis, the number 
and type of nuclei, and the precise muon site. With the 
orientation of P(0) fixed in the experiment, we used the 
approximation method to perform an exhaustive search 
for the function P(t) that best describes the ZF-/iSR sig- 
nal of La2- a; Sr a ;Cu04 using numerous combinations of 
the muon site and the surrounding nuclei. Figure Eta) 
shows P(t) calculated by the approximation method for 
a handful of the numerous potential muon sites consid- 
ered in our study. Some of these sites were considered 
in earlier fj.SK studiesj 34 ' 35 As demonstrated in Fig. OJb) , 
only nearest-neighbor nuclei need to be included in the 
calculation of P(t). In particular, note that both B a 
and £?b correspond to the same muon site, but the cal- 
culation of P(t) for B\y includes 3 additional Cu nuclei 
located further away from the muon (see Table fl} . The 
strong dependence of P(t) on the muon site means that 
ZF-/iSR can be used to accurately determine the location 
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FIG. 6: (Color online) Time evolution of the muon spin polar- 
ization P(t) calculated for the muon D site of Fig. [4] and dif- 
ferent orientations of the initial muon spin polarization P(0). 
The upper dashed curve shows P(t) for P(0) parallel to the 
c-axis. The lower solid curves show P(t) for different orienta- 
tions of P(0) in the a-b plane, labeled by the angle between 
P(0) and the a-axis. 



TABLE I: Muon sites of the polarization functions plotted in 
Fig. [5] The spatial coordinates of the muon site are denoted 
by multiples of the lattice constants a, b, and c of La2Cu04. 
The nearest-neighbor nuclei used in the calculation of P(t) 
for each muon site are listed in the third column. Two Cu 
nuclei are denoted as 2Cu, a La nucleus at site 1 (see Fig. [4]) 
is denoted by La(l), etc. 



Label 


Muon Site 


Nuclei 


M 


(0.5, 0.0, 0.096) 


2Cu, 2La(2) 


TI 


(0.2, 0.0, 0.15) 


2Cu, La(l), 2La(2) 


T2 


(0.225, 0.0, 0.225) 


2Cu, La(l), 2La(2) 


Ul 


(0.12, 0.0, 0.11) 


2Cu, La(l), 2La(2) 


A 


(0.0, 0.0, 0.212) 


5Cu, La(l), 4La(2) 


B a 


(0.1, 0.0, 0.1) 


Cu, La(l), 2La(2) 


B b 


(0.1, 0.0, 0.1) 


4Cu, La(l), 2La(2) 


D 


(0.12, 0.0, 0.219) 


Cu, La(l), 4La(2) 
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FIG. 7: (Color online) Comparisons between P(t) calcu- 
lated for the muon A and D sites and (a) the ZF-/^SR sig- 
nal of Lai.834Sro.i66Cu04 at T = 158 K measured with P(0) 
perpendicular to the c-axis, and (b) the ZF-/iSR signals of 
Lai. S7 Sro.i3Cu0 4 at T = 144 K from Ref. [3, for P(0) per- 
pendicular to the c-axis, and P(0) parallel to the c-axis. Note 
that the calculation for P(0) 1c assumes that the direction 
of P(0) makes an angle of 45° with the a-axis. Also, the 
ZF-/iSR signals shown here have been divided by the initial 
asymmetry a. 



of the implanted /i + in situations where there is a single 
muon site, the muon does not diffuse, and there are no 
additional sources of relaxation {e.g. electronic moments 
and/or loop-current order). 

Of all the muon sites we considered, only the polariza- 
tion function P(t) calculated for the site (0.120, 0, 0.219) 
shown in Fig. [4] (which we refer to here as the D site) 
accounts for the observed ZF-/J.SR spectra (see Fig. [7]). 
This site is located approximately 0.7 A from the api- 
cal oxygen, which agrees with the widely held view that 
the fi + bonds to an oxygen atom in cuprates. While it 
is not exactly one of the muon sites suggested in earlier 
workS ) 34 i 35 it is consistent with ZF-/1SR measurements of 
the antiferromagnetic phase of La2Cu04 by Hitti et al^ 



that restrict the muon site to the a-c plane. 

As shown in Fig. 0] and indicated in Table HI five La nu- 
clei of spin 1 — 7/2 and one Cu nucleus of spin 1 = 3/2 are 
included in the calculation of P(t) for the muon D site. 
Figure [6] shows how P(t) is dependent on the orientation 
of P(0) with respect to the crystal lattice. Although the 
orientation of P(0) in the a-b plane was random in our ex- 
periments, the ZF-/iSR spectra were reproducible when 
the sample was rotated about the c-axis to a different 
position. We can understand this as being a consequence 
of the equivalency of the sites (0.120, 0, 0.219) and (0, 
0.120, 219) in the tetragonal unit cell. Since these are oc- 
cupied by the muon with equal probability, the ZF-/xSR 
spectrum is an average of P(t) for these two sites. For 
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the case P(0) _L c, at any orientation this superposition 
is equivalent to P(t) calculated with P(0) lying in the 
a-b plane and making an angle of 45° with the a-axis. 
As shown in Fig.[7Ja), there is fairly good agreement be- 
tween the calculation assuming the muon is located at the 
D site and the ZF-/xSR spectrum of Lai.834Sro.i66Cu04 
measured with P(0) _L c. This is also shown to be the 
case for the A site, where the muon is positioned directly 
above the apical O. To distinguish between these two sites 
we consider the dependence on the angle between P(0) 
and the c-axis. In Fig. [Tfb) the calculations of P(t) for 
both P(0) _Lc and P(0) || c are compared to ZF-/1SR spec- 
tra of Lai. 8 7Sr .i3CuO4 from Ref. YL4. While this data is 
less accurate at late times due to lower muon counts, it is 
clear that the calculation for the A site fails to describe 
the ZF-/xSR signal with P(0) || c. 

IV. DISCUSSION AND CONCLUSIONS 

The ZF-/1SR measurements presented here and in 
Ref. [HI show no evidence for any form of static mag- 
netism that can be directly linked to the pseudogap phase 
of La2- ;c Sr ;c Cu04. Our determinaton of the muon stop- 
ping site enables a good estimate of the magnitude of 
the average local magnetic field that should have been 
detected if static orbital-like magnetic order of the kind 
observed by spin-polarized neutron diffraction"^ were 
present. The neutron results are in general agreement 
with the loop-current pattern of the O77 state proposed 
in Refs. |J and l36l. which consists of two oppositely circu- 
lating current loops per unit cell. If the orbital currents 
flow in the CuC>2 layers along the Cu-0 and 0-0 bonds as 
originally proposed, a [i + residing at the D site (0.120, 0, 
0.219) will experience a local field of 308 G//is. For the 
maximum possible ordered magnetic moment of 0.02 [Lb 
deduced from the spin-polarized neutron measurements 
of Lai.gi5Sro.o85Cu04fi^ the corresponding dipolar mag- 
netic field sensed by the muon is 6.2 G. In this case the 
damped ZF-//SR signal should oscillate with a period of 
~12 /is, such that one nearly complete oscillation is ob- 
served over the time range of our measurements. 

The moments of the orbital-like magnetic order ob- 
served in hole-doped cuprates by spin-polarized neutron 
diffraction are actually pointing at an angle of roughly 
45° with respect to the c-axis. This is compatible with 
the orbital currents of the O/j phase flowing out and 
back into the Cu02 layers through the apical oxygen 
atoms i2I In this scenario, the corresponding moments 
that are perpendicular to the faces of the CuOg octa- 
hedra in La2_ a; Sr 2: Cu04 point along directions making 
an angle of 61° away from the c-axis. For this arrange- 
ment with an ordered moment of 0.02 the average 
local field that a muon at the D site in La2-zSr :l: Cu04 
would detect is about 90 G. Since 0.02 /j,b is the maxi- 
mum possible value of the ordered moment observed in 
Lai.gi5Sro.o85Cu04, and the orbital-like ordered moment 
in other cuprates^— generally decreases with increased 



hole-doping, one might think of 90 G as an upper limit 
for the average local field in the x > 0.13 samples consid- 
ered here. On the other hand, a larger value is possible 
if the ordered moment in Lai.9i5Sro.o8sCu04 is reduced 
due to a higher degree of disorder and/or competition of 
the orbital-like magnetic order with the spin-density wave 
(SDW) order present in lower doped samples. Competi- 
tion with SDW order has been suggested as a potential 
explanation for the smaller ordered moment observed in 
YBa 2 Cu3 06.45 compared to that of YBa 2 Cu306.5i£ Re- 
gardless, the contribution of orbital-like magnetic order 
to the ZF-/xSR signal will be strongly damped if it is 
short range as in the case of Lai.gisSro.ossCuO.!. Yet no 
such component is observed in any of the La2- a; Sr :I ;Cu04 
samples we measured. 

The failure here and in Ref. [H] to detect orbital-like 
magnetic order in La2^ 2; Sr 2 ;Cu04 of the kind observed by 
spin-polarized neutron diffraction may indicate that the 
local fields are rapidly fluctuating oustide the /iSR time 
window. It has also been suggested that the may de- 
stroy loop-current order in cuprates, 38 although thus far 
it has also eluded detection by NMR/NQR. On the other 
hand, since the polarized neutron diffraction experiments 
cannot deduce magnetic volume fractions, the orbital-like 
magnetic order could still be associated with a small mi- 
nority phase that evolves with hole doping. 13 Having said 
all of this, there is currently no disagreement between the 
polarized neutron diffraction and ZF-/iSR experiments on 
La2- 2; Sr 2 ;Cu04, since there is no overlap in the doping 
range of the samples studied by these two techniques. 
The ZF-/iSR experiments have focussed on x > 0.13 sam- 
ples to avoid major contributions from static SDW or- 
der or spin-glass-like magnetism. Given that the orbital- 
like magnetic order observed by spin-polarized neutron 
diffraction in Lai.gi5Sro.os5Cu04 is quite weak and has 
not been observed in other cuprates beyond p = 0.135, 
it is conceivable that it is not present in x > 0.13 sam- 
ples. What can be said is that the ZF-/J.SR measurements 
of La2-a;Sr x Cu04 in the strontium concentration (hole- 
doping) range 0.13 < x < 0.19 do not support theoreti- 
cally predicted loop-current phases, and hence favour an 
alternative explanation for the unusual magnetic order 
detected by spin-polarized neutron diffraction at lower 
hole doping. 
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